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SUMMARY 

We analyse the linear thermal stability of a 2D periodic structure (alterna- 
tively hot and cold) in a uniform magnetic field. The energy equation in- 
cludes wave heating (assumed proportional to density), radiative cooling 
and both conduction parallel and orthogonal to magnetic lines. The equili- 
brium is perturbed at constant gas pressure. With parallel conduction 
only, it is found to be unstable when the length scale 1// is greater than 
45 Mm. In that case, orthogonal conduction becomes important and stabilizes 
the structure when the length scale Ij^ is smaller than 5 km. On the other 
hand, when 1J^ is greater than 5 km, the thermal equilibrium is unstable, 
and the corresponding time scale is about 10^ s : this result may be 
compared to observations showing that the lifetime of the fine structure 
of solar prominences is about one hour ; consequently, our computations 
suggest that the size of the unresolved threads could be of the order of 
10 km only. 

FUNDAMENTAL EQUATIONS OF THERMAL EQUILIBRIUM. 

We use the 2D "chessboard" of figure 1 , which displays a periodic hot 
(Tl) and cold (T2) structure in a uniform magnetic field (this could be 
the case in the central parts of prominences, see Leroy et al., 1983). 

We write equations for thermal equilibrium of hot and cold cells, as : 



are respectively the heat flux parallel and orthogonal to the magnetic 
field, pi and p 2 are the densities of respectively hot and cold cells. 
k Q // and k 0 ^ are conduction coefficients ; k 0- Ldepends on the strength 
of the magnetic field B. 
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We assume that the gas pressure remains constant = f>2 ^2 ~ = constant 

and use the cooling function Q (T) given by Hildner (1974). Unknown quantities 
are h, Tj, T 2 . The equilibrium state (h, T 2 ), when Ti is fixed, is given by 
the set of 2 equations above. Possible solutions are shown in figure 2 (top). 

ANALYSIS OF THERMAL STABILITY 

We perturb the equilibrium at constant gas pressure P and heating h (T-^T+^T). 
Equations are linearized assuming that $ T e./*^ , where J& is a growth rate. 

We assume also that there is no motion in a direction perpendicular to the 
magnetic field, so that lj_ = constant. Hence, mass conservation gives !///> = 
constant. We get two solutions forjS ; the first one is always negative (stable 
solution), but the second one can either be positive (unstable) or negative 
(stable). It depends on the values of equilibrium parameters 1//, 1J, and T^. 

The magnetic field strength B was kept constant (1 Gauss). 

RESULTS 

Figure 2 (bottom) gives the growth rate as a function fo l^and hot tempera- 
ture T 1 (parallel conduction was neglected there). 

The thermal equilibrium is unstable when £> 0 : this is always the case when 
Ij. is too large (> 10 6 m). When lj^ is smaller than 10^ m, the equilibrium may 
be stable if T| does not exceed a critical value Tmax (Tmax = 10^ K for 
ljL = 10^ m) . 

Figure 3 gives^as a function of 1^ and 1// for a 2 D model (T]^ was kept 
constant). It shows that, when 1// < 45 Mm, the thermal equilibrium is stable ; 

when 1// > 45 Mm, it is unstable, unless IX < 5 km. When IX > 5 km, the time 

scale for instability is approximately equal to 10^ s and corresponds to the 
observed life time of the fine structure in solar prominences. This result 

suggests that the size of thin threads could be as small as 10 km. 

The temperature of hot cells (T^) used in the computations was 10^ K. 
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figure 1 : the geometry of the model : the structure is periodic in x direction 
(length scale lj^) and y direction (length scale 1//). The magnetic field B is 
parallel to y. z is the vertical axis. White areas are hot (T^) and tenuous 
( O i) ; dashed cells are cold (T 2 and dense (# 2 ^ • 
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top : log Q (T) + (orthogonal conduction) / ^ as a function of log (T) for 
differnt values of the length scale lj^ ( ID calculation with lj. = 10®, 

10 ^ 5 10 5 , 10 * and 10* m). Equilibrium solutions are located at points 
A, B, C (intersection with the straight line log (T) + log (AJi )). 

When the heating becomes too large, the cold solution A does not exist 
any more. 


bottom : the growth rate as a function of hot temperature (1 D calculation) 
The function sgn (ft) / log |ftl is displayed for different values of lj_ 
(same as above). J 
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figure 3 : the growth rate for the 2 D calculation as a function of 1// and 
IX • Isocontours of sgn ( log IfM are displayed. The dotted area is stable 
( jS> < 0) ; the white one (top) is unstable ( > 0). The dashed region (left 
and bottom) represents the domain where a cold equilibrium does not exist 
(see figure 2). 
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